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ABSTRACT 

We combine photometric observations of high-redshift [z > 4) quasars, obtained at 
submilhmetre (submm) to miUimetre wavelengths, to obtain a mean far-infrared (FIR) 
(rest- frame) spectral energy distribution (SED) of the thermal emission from dust, 
parameterised by a single temperature (T) and power-law emissivity index (/?) . Best- 
fit values are T = 41 ± 5K, /3 — 1.95 ± 0.3. Our method exploits the redshift spread 
of this set of quasars, which allows us to sample the SED at a larger number of 
rest wavelengths than is possible for a single object: the wavelength range extends 
down to ~ 60 fj,m, and therefore samples the turnover in the greybody curve for these 
temperatures. This parameterisation is of use to any studies that extrapolate from a 
flux at a single wavelength, for example to infer dust masses and FIR luminosities. 

We interpret the cool, submm component as arising from dust heated by star- 
formation in the quasar's host galaxy, although we do not exclude the presence of dust 
heated directly by the AGN. Applying the mean SED to the data, we derive consistent 
star-formation rates ~1000Mq yr~^ and dust masses ^ 10^ Mq, and investigate a 
simple scheme of AGN and host-galaxy coevolution to account for these quantities. 
The timescale for formation of the host galaxy is 0.5-lGyr, and the luminous quasar 
phase occurs towards the end of this period, just before the reservoir of cold gas is 
depleted. Given the youth of the Universe at z = 4 (1.6Gyr), the coexistence of a 
massive black hole and a luminous starburst at high redshifts is a powerful constraint 
on models of quasar host-galaxy formation. 
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1 INTRODUCTION the dust obscuration necessarily makes the sources optically 

faint. 



Cosmological studies at FIR-submm wavelengths are an es- 
sential complement to optical investigations, in that the dust 
that extinguishes the optical/UV light from AGN or mas- 
sive stars is that which thereby shines in the infrared. This 
effect should not be underestimated: the integrated energy 
density in each of the FIR and optical extragalactic back- 
ground light is comparable (Paget et al. 1996; Hauser et 
al. 1998). It has been suggested that cosmological submm 
sources are high-redshift equivalents of ULIRGs, their lumi- 
nosity deriving from bursts of dust-obscured star-formation 
accompanying the assembly of massive elliptical galaxies, in 
the course either of gas-rich hierarchical major-mergers or ab 
initio monolithic collapse (Barger, Cowie & Sanders, 1999, 
and references therein). However, such speculations are dif- 
ficult to confirm via direct optical follow-up, because of the 
large beam-sizes of current submm telescopes, and because 



An alternative approach is to study the submm proper- 
ties of well-characterised classes of objects, whose precisely- 
known positions and redshifts, physical properties and pop- 
ulation statistics, facilitate follow-up observations and inter- 
pretation. To this end, we have been studying dust contin- 
uum and CO emission from luminous {Mb ^ —27.5), high- 
redshift {z J> 4) quasars. At low redshifts, the hosts of lumi- 
nous quasars are invariably massive elliptical galaxies (Bah- 
call et al. , 1997; McLure et al. , 1999). In the local Universe, 
massive black holes, fossil remnants of AGN now starved of 
fuel, are found in the cores of most galactic bulges; more- 
over, their mass is correlated with the stellar mass of the sur- 
rounding spheroid, Mbh = 6 x 10~^A/sph (Magorrian et al. , 
1998). This correlation hints at a deep connection between 
black-hole build-up through accretion and star-formation in 
the host galaxy. Attributing the submm-luminosity of high- 
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Table 1. Summary of SED parameters used by other authors 



Authors 


Sample 


ifir 


T 


/3 


Dunne et al. 


SLUGS 


10" 


36 


1.3 


Carico et al. 


LIRGs 




32 


1.6 


Lisenfeld et al. 


LIRGs 


IQll.S 


38 


1.6 


Benford et al. 


High-z objects 


1013 


52 


1.5t 


Hughes et al. 


High-2 objects 


1013 


50t 


1.5t 


McMahon et al. 


2 > 4 quasars 


1013 


50t 


1.5t 


This work 


z > A quasars 


10" 


42 


1.9 



t values presumed, not determined 



redshift quasars to star-formation therefore accords with the 
speculations about the origins of the submm sources. 

Before one can estimate physical quantities, such as 
dust temperature and star-formation rate, it is necessary 
to characterise the form of the dust-emission spectrum of 
these objects. It is also important to quantify variations of 
the SED with object class, luminosity or redshift, other- 
wise, if one adopts a single template for all sources, it is 
not legitimate to adopt more than one derived parameter 
as an indicator of physical status. In fitting thermal SEDs, 
however, we are constrained by the paucity of photometric 
data, and often it is inappropriate to fit both the emissivity 
and the temperature as free parameters, when we have only 
three or four photometric points per source. Previous work, 
which considered object-by-object temperature fits, adopted 
an assumed value for the emissivity index (Bufi^ey et al. , in 
prep.). In contrast, the purposes of this Letter are, first, to 
investigate the possibility of co-adding data from a sample 
of luminous, 2^4 quasars, to find a self-consistent, over- 
all best-fit temperature and emissivity; and to characterise 
the uncertainties involved in doing so, whether they are due 
to experimental errors or to physical differences. Finally, we 
discuss applications of the mean SED, and try to justify the 
inferred parameters in terms of a simple evolutionary model. 

Throughout, we assume a A-dominated cosmology 
nM=0.3, nA=0.7, ifo=65kms"^ Mpc"\ 



2 MEAN ISOTHERMAL SPECTRUM 

2.1 Data 

The 2 > 4 datasets are from our work with IRAM (1.25 mm) 
and JCMT (SCUBA 850 fim, 750 ^m, 450 ^m and 350 ^m; 
UKT14 800 jj,m) as described in companion papers (Isaak et 
al. 1994; Omont et al. , 1996a; Buffey et al. , in prep.); and 
350 pm CSO photometry is from Benford et al. (1999). Ad- 
ditionally, we include the z = 3.91 quasar APM08279-I-5255 
(Lewis et al. 1998). All the sources have been detected at 
least four different wavelengths. An important caveat is that 
the quoted statistical uncertainties do not reflect the sys- 
tematic errors due to flux calibration. Therefore, we assume 
calibration errors, combined in quadrature with the random 
uncertainties, of 20 percent. 

2.2 Assumptions 

The rest-frame far-infrared (FIR) SED can be described by 
a modified blackbody spectrum, assuming that the grains 
are in equilibrium with the radiation field and that the dust 



cloud is optically-thin at submm wavelengths; S„ oc — . 

efcT -1 

Note, however, that some ultraluminous starbursts (such 
as Arp220; Rowan- Robinson & Efstathiou, 1993) are opti- 
cally thick at the shortest rest-wavelengths sampled by the 
present data (60/im). This method determines the most ef- 
fective isothermal parameterisation of the long-wavelength 
emission, a simple, convenient description that is frequently 
adopted {e.g. Benford et al., 1999; Lisenfeld, Isaak & Hills, 
2000; Dunne et al., 2000), though in reality, a distribution 
of dust temperatures is likely. Shorter wavelength radiation, 
in the near-to-mid infrared region (~10/im rest-frame), is 
likely to be contaminated by the hotter dust that absorbs 
the harder continuum due to the AGN itself {e.g. Rowan- 
Robinson, 1995). The contribution from star formation, on 
the other hand, emerges in the submillimetre region. To as- 
sess the contamination of the submillimetre luminosity due 
to the AGN component requires better mid-infrared con- 
straints than exist for all but a few high-redshift quasars 
(see, for example, Rowan-Robinson, 2000). We therefore 
simply assume that the component extracted by our flt- 
ting procedure is appropriate for the extrapolation of a far- 
infared luminosity (Lfir) from which a star formation rate is 
to be calculated. We caution, however, that the unforeseen 
presence of a hot component could skew the fit to higher 
temperatures, leading us to overestimate the star formation 
rate. 



2.3 Method 

Adopting T and /3 as free parameters, we employ a 
X^-minimisation procedure to determine the best single- 
component fit to the data, which is coadded in a self- 
consistent way. A subtlety arises because for a given ob- 
served wavelength, the range of objects' redshifts causes us 
to sample their spectra at different rest wavelengths. There- 
fore, if we choose to normalise the data at a given rest wave- 
length, a spectrum must be assumed so that an extrapola- 
tion to this reference from the nearest observed data point 
can be performed. The best-fit spectrum thus determined 
then itself becomes the template for the extrapolation, and 
the process is iterated. In practice, this procedure converges 
to a solution rapidly {i.e. without roaming around parameter 
space) so for these datasets at least, we conclude that the 
overall x^-minimised fit is, indeed, the best self-consistent 
solution. 



2.4 Results 

The contours are shown in Figure 1, while the renor- 
malised data and best-fit SED axe plotted in Figure 2. 
The best self-consistent fit corresponds to /? = 1.95 ± 0.3, 
T = 41 ± 5K (Icr errors). A similar procedure applied to 
the low-redshift, less-lumimous samples of LIRGs (Lisen- 
feld, Isaak & Hills, 1999) and local galaxies (Dunne et al. , 
2000) gives results inconsistent at the > 2a level (Figure 1). 

We can use the contours to predict the outcome of 
curve-fitting in which 13 is held fixed, while the temperature 
is varied. For other studies of high-redshift quasar SEDs, 
the /3 = 1.5 adopted by Benford et al. (1999) and Buffey 
et al. (in prep.) would give r=50K as mean best-fit, which 



© 0000 RAS, MNRAS 000, 000-000 



The submm spectral energy distribution of high-redshift quasars 3 



3.0 



2.0 



1.0 



10.0 F 




z>4 quasars 
LIRGs 

Local galaxies 



40 60 
Temperature, K 



Figure 1. Contours of x^i in /3 — T space, for our combined 
photometry of z > 4 quasars. The star-shaped symbol indicates 
the best-fit position, and both 1— and 2-a contours are plotted 
(solid). For comparison, the best-fit and 1- and 2-ct contours for 
local LIRGs (from Lisenfeld, Isaak & Hills, 1999) (asterisk and 
dashed), and for local galaxies (from Dunne et al. , 2000) (cross 
and dotted), are also shown. The dot-dashed line is a locus of con- 
stant luminosity derived from the best-fit T and /3, for constant 
850 /im flux at z = 4.5. 



is indeed close to that found by those authors, and consis- 
tent with the T=50K assumed by McMahon et al. (1999). 
Variation in luminosities (hence, eventually, star-formation 
rates) estimated from such best-fits do not vary substan- 
tially, because Lfir cx (3 -I-/3)! x T^+i^, whilst the best-fit T 
is monotonically decreasing with (3: the dot-dashed line in 
Figure 1, a locus of constant luminosity, illustrates this. This 
is to be oxpoctod, since luminosity is an integral under the 
curve defined by the data points. On the other hand, using 
best fits to the datasets consisting of other objects (listed in 
Table 1), would give a discrepancy in luminosity of up to a 
factor three at the redshifts of our quasars. 



3 APPLICATIONS 

In this section, we discuss some applications for observation, 
and implications for interpretation, of this mean SED. 



3.1 Millimetre— submillimetre spectral index 

The stoop slope of cool, thermal SEDs over the mm-submm 
range makes knowledge of the precise spectral shape impor- 
tant, when comparing data obtained with different instru- 
ments operating at different wavelengths. Consider, for ex- 
ample, the spectral index between 850 /xm and 1.25mm, ap- 
propriate for converting between ,]CMT/SCUBA and IRAM 
data, respectively. The ratio of fluxes as a function of red- 
shift for our best fit to «> 4 quasars, together with the fits 
to LIRGs and local galaxies, is shown in Figure 3. For ease 



1.0 r 



0.1 r 




100 


: 

















100 

Rest wavelength, /xm 



Figure 2. The best-fit, rest-frame SED, with renor- 
malised data points. Each separate symbol denotes a dif- 
ferent quasar: BR1202-0725 (filled circles), APM08279-t-5255 
(open circles), BRI1335-0417 (squares), BR1117-1329 (stars), 
BRI0952-0115 (upright triangles), BR1033-0327 (inverted tri- 
angles), BRl 144— 0723 (diamonds). The solid line corresponds to 
T = 41K, /3 = 1.95, the x^-niininiiscd best fit to these points, 
while the dashed and dotted lines represent the best fits to the 
data of Lisenfeld ct al. and Dunne ct al. (as in Figure 2). The 
dot-dashed line is the starburst model of Rowan- Robinson & Ef- 
stathiou (1993). The curves are normalised at 850/im/(l-|-4.5). 




Redshifl 

Figure 3. Ratio of observed flux at 850 ;um to that at 1250 /im — 
i.e. , the conversion between typical observations made with 
JCMT/SCUBA and with IRAM respectively. The solid lines are 
best fit (centre) and ±1- and 2— a . The apparent convergence at 
2 ~ 6 is due to the tight bunch of data points near rcst-framc 
150 ^m (see Figure 2). The dashed and dotted lines correspond, 
as in Figures 1 and 2, to the best- fits for the Lisenfeld et al. and 
Dunne et al. samples, whilst the dot-dashed line depicts Arp220. 
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Table 2. 2 > 4 sources used in the fit, from Omont et al. (1996), with measured and extrapolated properties based upon the 850 /im 
fluxes of Buffey et al. (in prep.) 



Object 


z 


Mb 


10^ M0 


1013 L0 


10^ M0 


10" M0 


10" M0 


SFR 

103 Ma 
yr 


10^ yr 


2^max 




BRI0952-0115 


4.43 


-28.0 


5.5 


2.0 


1.7 


1.4 


9.1 


1.2 


0.9 


9.5 


2.9 


BR1202-0725 


4.69 


-28.8 


11 


5.8 


5.0 


4.1 


19 


3.5 


0.7 


7.9 


3.3 


BR1335-0417 


4.40 


-27.6 


3.8 


2.0 


1.7 


4.3 


6.3 


1.2 


0.6 


6.9 


3.2 


BR1033-0327 


4.50 


-27.9 


5.0 


1.0 


0.8 


<0.9t 


8.3 


0.6 


1.7 


00 


2.3 


BR1117-1329 


3.96 


-28.2 


6.6 


1.9 


1.7 


<1.2t 


11 


1.1 


1.2 


10 


2.5 


BR1144-0723 


4.14 


-27.7 


4.1 


1.0 


0.9 


<1.3t 


6.9 


0.6 


1.4 


20 


2.4 


Median Quasar^ 


4.5 


-28.0 


5 


2 


2 


1.5 


8 


1.2 


0.7 


7.4 


3.2 



Note: All derivations assume Qm = 0.3, Qa = 0.7, Ho = 65kms^l Mpc~l 
t Scr upper limits 

t represents median of the observed quasars, consistent with an object having Sgso ~ 15mjy and Mg = —28.0. 



of reference, a linear approximation to the ratio for our best 
fit, good out to z - 7, is R{-^) « 3.7 - z/A. 



3.2 Physical characteristics 

Table 2 lists physical properties inferred for each of the 
quasars contributing to the fit. Dust masses (Md) and far- 
infrared luminosities (Lfir) are calculated as in McMahon 
et al. (1999), though with our new T and /3. The final stellar 
mass of the host spheroid (Mgph) is estimated using the lo- 
cal relation between bulge and black hole mass (Magorrian 
et al. , 1998), with the latter calculated from the absolute 
i?-band magnitude (Mb) via simple assumptions about ac- 
cretion (Section 3.3), and a bolometric correction of 12 after 
Elvis ct al. (1994). Molecular gas masses (Mhj) are calcu- 
lated from observed CO emission-line luminosities (Omont 
et al. 1996b; Guilloteau et al. 1997, 1998), using a conver- 
sion factor Mna/ico = 4.5M0(Kkms"i pc^)"^ (Solomon, 
Downes & Radford, 1992); though we caution tliat a cali- 
bration as much as five times smaller than this has been 
reported to be more appropriate for some ULIRGs (Downes 
& Solomon, 1998). 

The substantial star-formation rates and far-infrared 
luminosities might cause us some concern: Lfir ~ Lb, so 
we should consider whether the AGN itself is an impor- 
tant contributor to the cool dust luminosity. Of course, 
the equally-substantial quantities of of cool molecular gas 
(Mh2 > Mg ) make ideal nests for hatcliing young stars; 
and stars are required to manufacture the large mass of 
dust that is observed (Md > 10^ M©). Moreover, the 00(5- 
4)/1.35mm-continuum of BR1202— 0725 has been resolved 
into two components (Omont ct al. , 1996), only one of 
which coincides with the quasar, the other showing no sign 
of AGN activity. Yun et al. (2000) have detected the three 
CO-luminous quasars at 1.4 GHz: the reported fluxes are 
consistent with synchrotron emission from supernova rem- 
nants, and imply star-formation rates ^lOOOMgyr^i, in 
agreement with the submm predictions (Section 3.3, and 
Table 2). Therefore, although we do not exclude the possi- 
bility that some of the submm flux is due to AGN-heating, 
here we shall explore the hypothesis that it derives wholly 
from star-formation. 

This leaves us with somewhat extraordinary objects, 
which are not only undergoing intense bursts of dust- 



obscured star formation, but simultaneously liosting power- 
ful quasars. The age of the Universe at 2; = 4 is 1.6 Gyr: this 
lack of time imposes severe constraints upon the processes 
governing both central black-hole growth (Turner 1993) and 
the chemical evolution necessary to supply the large ob- 
served dust masses (Edmunds & Eales 1998). The statis- 
tics are, currently, very poor (only three z > 4 objects with 
CO-detections, biased towards the most luminous objects), 
nevertheless we assume that the Median Quasar of Table 
2 epitomises its class, to guide our quest for a consistent 
hypothesis. 



3.3 Quasar and host galaxy coevolution 

The far-infrared luminosity (Lfir) can be used to esti- 
mate the current star-formation rate, given a stellar Ini- 
tial Mass Function (IMF) and a prior star-formation his- 
tory: 5Fi?/ Moyr^i = 'I'Lfir/IO"' Lg. To estimate *, we 
use the stellar population synthesis model of Leitherer et 
al. (1999), STARBURST99, assumlug a constant star-formation 
rate and that Lfir = Lboi. A significant cause of uncertainty 
is due to the choice of IMF: introducing faint, low-mass 
stars has the effect of increasing the total mass while hardly 
affecting the luminosity, therefore one must infer a higher 
star- formation rate for a given Lboi. We adopt a Salpeter 
IMF {4>{m) oc m~^'35^ Salpeter 1955) within the mass range 
1 Mq < m < 100 M0 . ^' depends much more weakly upon 
metallicity, and upon tlic age of tfie burst for lifetimes longer 
than ~ lO'^ yr. We take = 0.6, which corresponds to a 
metallicity Z = Zq and an age t « 0.5 Gyr. Another source 
of uncertainty is the fraction of starlight that is intercepted 
by dust; liere we make the simplifying assumption that this 
fraction is unity. Star-formation rates thus derived are pre- 
sented in Table 2. 

The spheroid-formation timescales in Table 2 corre- 
spond to the total time required to convert the inferred 
spheroid masses from gas into stars, at constant star- 
formation rate; not forgetting gaseous material expelled, 
from massive stars, as winds and supernovae (we assume 
a return fraction of ~20 percent). It is through this outflow 
process tiiat the interstellar medium (ISM) becomes chemi- 
cally enriched: if half the mass in metals condenses into dust, 
then chemical yields derived from stellar evolution models 
[e.g. Maeder, 1992) predict that ~1 percent of the mass pro- 
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cessed through stars (for our IMF) is returned to the ISM 
in the form of dust — neglecting however the intricacies of 
dust-destruction and -dispersion by supernovac. Thus, af- 
ter 0.5 Gyr, the stellar mass of the Median Quasar would 
be ^6x10^^ Mq, its gas mass is ^2x10"'^^ Mq, and its dust 
mass is «3xlO®M0 — near the maximum it can attain for 
these parameters, before becoming recycled into long-lived 
stellar remnants (Edmunds & Bales, 1998). 

Meanwhile, to describe the evolution of the AGN, we 
adopt the simplest assumptions about accretion onto a mas- 
sive black hole (Rees, 1984): that the mass flow Mace is reg- 
ulated by radiation pressure — i.e. accretion occurs at the 
Eddington rate — and that the accretion efficiency e ~ 0.1, 
so 

Mbh 



Lhoi = eMa, 

where ctt = 
and rup = 
rate is Ma 



ATrGrUpCMbh 



3 X 10 



13 



109 Mo 



J©, 



(1) 



1.67 X 10 

_ _22 

<= - e/0.1 109 M(7 



CTT 
25 , 

g is the proton mass. The accretion 
Mbh^ ]y[g yr~^ and the mass e-folding 



o 



eAf, 



bh<: 



time for the growth of the black hole is Tbh ~ 
^ X 4.5 X 10*^ yr. Supposing that Eddington-limited accre- 
tion starts to apply above a seed black hole mass of about 
10® Mq (Haehnelt, Natarajan & Rees, 1998), it would take 
~0.5 Gyr to build the required quasars — close to the star- 
forming lifetime of the host galaxy, discussed above. 

The total (baryonic) mass of the galaxy has been 
inferred from the black hole mass by taking the local 
remnant black hole-bulge correlation (Magorrian et al. . 
1998). Critics may object that we do not measure the fi- 
nal mass of the black hole, because we observe it while it 
is still actively accreting. However, by the time the black 
hole's mass is 5 x 10^ Mq , its Eddington accretion rate 
is lOOMoyr^^and rising, and it is competing with star- 
formation (~1000 M0 yr^^ ) for the rapidly-diminishing sup- 
ply of cold gas. In the absence of inflow, this cannot be sus- 
tained for long. Therefore, we obtain one consistent scheme 
by assuming that the black holes' masses do not grow much 
greater than the values deduced from the quasar light. This 
is in accord with the measurements of remnant black hole 
masses in the local Universe, which give at most 10®'® M©. 

The quantities ^imax and Zmin in Table 2 bracket the red- 
shift range over which the quasars are forming stars, based 
upon the inferred spheroid mass and constant star-formation 
rate. The severity of the time constraint at high redshifts 
is dramatically illustrated, especially for BR1033— 0327, 
which has a star-forming lifetime longer than the age of 
the Universe at its rcdshift. In this case, we could re- 
solve the dilemma by invoking a different IMF and star- 
formation/accretion histories: for example, exponential gas- 
consumption and a phase of super-Eddington accretion, or 
more low-mass stars leading to a higher star-formation rate, 
would give rather shorter lifetimes. We refrain from investi- 
gating this here, however, for lack of space and datar volume, 
but we shall examine the problem in later papers. 



4 CONCLUSIONS 

We have constructed a mean SED, over the rest-frame wave- 
length range 60-300 fj,m, for luminous z > 4 quasars, assum- 
ing an isothermal dust spectrum. We consider applications 



of this SED, such as the conversion between commonly- 
used submm/mm wavebands, and the estimation of star- 
formation rates and dust masses. Using the observed local 
remnant black hole to bulge mass relation to estimate the 
initial gas mass, we deduce star-formation timescales and 
sketch a simple evolutionary scenario. 

Estimation of star-formation rate and dust and gas 
masses obtained through submm-mm observations are po- 
tentially important constraints on models of quasar-host- 
galaxy cocvolution. However, the correct interpretation of 
such information requires detailed knowledge of the SED 
and its variation with parameters such as redshift and lumi- 
nosity. Otherwise, it is misleading to treat, as independent, 
quantities whose derivation follows from a single data point. 
Although we have here determined a "mean" SED, we do 
not suggest that all high-redshift quasars possess a unique 
dust temperature: much of the scatter about the mean is 
very likely real. The data are not yet extensive enough to 
justify a more detailed treatment, particularly since only 
the most luminous objects are suitable for multi-wavelength 
photometry. 

We emphasise, too, that it is important to justify the 
large derived star-formation rates and dust masses in terms 
of the evolution of stellar populations and chemical abun- 
dances: the simple treatments frequently adopted are not 
necessarily misleading, but neither are they immune to jus- 
tified criticism. Employing a skeletal evolutionary scheme, a 
consistent account emerges, within the many uncertainties, 
if we identify submm-luminous, high-redshift quasars with 
an early, gas-rich phase in the formation of massive galax- 
ies, at about the time when the dust mass reaches its max- 
imum, after the AGN has grown sufficiently luminous for 
it to be identified as a bright quasar and just before star- 
formation hcis depleted the fuel-supply upon which it sub- 
sists. We estimate that this quasar-starburst phase is short, 
lasting for little longer than 0.5 Gyr, though by changing 
our assumptions — for example, regarding the IMF or early 
accretion and star-formation history — , further consistent 
models could be constructed. Forthcoming observational re- 
sults (Isaak et al. , in prep.) will improve statistics and per- 
mit refinement of these speculations. 



ACKNOWLEDGMENTS 

RSP thanks PPARC for support, and RGM thanks the 
Royal Society. We thank the referee for comments. 



REFERENCES 

Bahcafl J.N., Kirhakos S., Saxe D.H., Schneider D.P., 1997, ApJ, 
479, 642 

Barger A.J., Cowie L.L., Sanders D.B., 1999, ApJ, 518, L5 

Downes D. & Solomon P.M., 1998, ApJ, 507, 615 

Dunne L., Ealcs S., Edmunds M., Ivison R., Alexander P., 

Clements D.L., 2000, MNRAS, 315, 115 
Edmunds M.G, Eales S.A., 1998, MNRAS 299, L29 
Elvis M. ct al., 1994, ApJS, 95, 1 

Guilloteau S., Omont A., McMahon R.G., Cox P., Petitjean P., 

1997, A&A, 328, LI 
Guilloteau S., Omont A., Cox P., McMahon R.G., Petitjean P., 

1999, A&A, 349, 363 



© 0000 RAS, MNRAS 000, 000-000 



6 R.S. Priddey & R.G. McMahon 



Hachnclt M.G., Natarajan P., Rees M.J., 1998, MNRAS, 300, 817 
Hauscr M.G. ct al., 1998, ApJ, 508, 25 

Hughes D.H., Dunlop J.S., Rawlings S., 1997, MNRAS, 289, 766 
Isaak K.G., McMahon R.G., Hills R.E., Withington S., 1994, MN- 
RAS, 269, L28 
Leitherer C. et al. , 1999, ApJS, 123, 3 

Lewis G.F., Chapman S.C., Ibata R.A., Irwin M.J., Totten E.J., 

1998, ApJL, 505, 1 
Lisenfeld U., Isaak K.G., Hills R., 2000, MNRAS, 312, 433 
Maeder A., 1992, A&A, 264, 105 
Magorrian J. et al. , 1998, AJ, 115, 2285 

McLure R.J., Kukula M.J, Dunlop J.S., Baum S.A., O'Dea C.P., 

Hughes D.H., 1999, MNRAS, 308, 377 
McMahon R.G., Omont A., Bergeron J., Kreysa E., Haslam 

G.G.T., 1994, MNRAS, 267, L9 
McMahon R.G., Priddey R.S., Omont A., Snellen I., Withington 

S., 1999, MNRAS 309, LI 
Omont A., McMahon R.G., Cox P., Kreysa E., Bergeron J., Pajot 

F., Storric-Lombardi L.J., 1996a, A&A, 315, 1 
Omont A., Pctitjcan P., Guillotcau S., McMahon R.G., Solomon 

P.M., Pccontal, E., 1996b, Nature, 382, 428 
Puget J.-L. ct al., 1996, A&A, 308, L5 
Rees M.J., 1984, ARA&A, 22, 471 
Rowan-Robinson M., 1995, MNRAS, 272, 737 
Rowan-Robinson M., 2000, MNRAS, 316, 885 
Rowan-Robinson M. & Efstathiou A., 1993, MNRAS, 263, 675 
Salpeter E.E., 1955, ApJ, 121, 161 

Yun M.S., Carilh C.L., Kawabe R., Tutui Y., Kohno K., Ohta K., 
2000, ApJ, 528, 171 



© 0000 RAS, MNRAS 000, 000-000 



